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The  purpose  of  this  3tudy  is  to  provide  an  analysis  of  hindcast  and  mea¬ 
sured  wave  and  wind  data  to  estimate  the  long-term  extreme  wave  heights  for 
the  St.  Paul  Harbor  area  at  Kodiak,  Alaska.  This  study  included  the  develop¬ 
ment  of  the  deepwater  extreme  wave  climate,  the  derivation  of  a  wave  energy 
attenuation  factor  for  deep-  to  shallow-water  transition  of  extreme  waves, 
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20.  ABSTRACT  (Continued): 

and  an  analysis  of  collateral  wind  data  to  support  the  results  of  the  analysis 
of  hindcast  wave  data. 

The  Wave  Information  Study  (WIS)  hindcast  data  applicable  to  the  Kodiak 
area  was  surveyed ‘for  extreme  wave  conditions.  This  survey  resulted  in  a 
sample  of  62  significant  wave  heights  for  the  years  1956  through  1975  (Appen¬ 
dix  A).  The  62  significant  wave  heights  were  scaled  and  plotted  according  to 
several  extreme  probability  models.  The  Extremal  Type  I  or  Fisher-Tippett 
Type  I  model  demonstrated  a  better  fit  with  the  data  than  did  the  other  models 
(Figure  1).  Extrapolated  significant  wave  heights  based  on  the  Extremal 
Type  I  model  are  listed  versus  return  period  in  Table  7.  The  extreme  wave 
predictions  computed  using  monthly  maximum  significant  wave  height  values  from 
measured  wave  data  (Table  9)  tended  to  be  smaller  than  the  hindcast  predic¬ 
tions.  This  difference  was  best  explained  by  the  fact  that  the  measured  data 
covered  a  time  span  of  about  2  years;  whereas,  the  hindcast  data  represented 
20  years  of  wave  climate.  The  longer  record  had  a  better  chance  of  capturing 
a  representative  sample  of  extreme  events  than  the  short  record. 

Extreme  wave  heights  computed  with  long-term  wind  data  provided  by  the 
National  Weather  Service  agreed  well  with  results  from  the  analysis  of  WIS 
hindcast  data  and  supported  the  overall  validity  of  the  study. 

Measured  wave  data  provided  by  the  Alaska  Coastal  Data  Collection  Pro¬ 
gram  were  analyzed  to  obtain  a  wave  energy  attenuation  factor  for  the  long- 
period  swell  crossing  the  reefs  into  St.  Paul  Harbor  at  Kodiak.  The  atten¬ 
uation  factor  represented  the  fraction  of  deepwater  wave  energy  that  is  not 
dissipated  in  the  transition  across  the  reefs  into  St.  Paul  Harbor.  The  at¬ 
tenuation  factor  had  a  value  of  0.374;  therefore,  a  significant  wave  height 
of  1  m  in  deep  water  corresponded  to  0.374  m  inside  the  harbor. 

The  long-term  extreme  wave  height  prediction  from  the  analysis  of  hind¬ 
cast  data  was  adjusted  to  represent  the  long-term  wave  climate  inside  St.  Paul 
Harbor  by  means  of  the  attenuation  factor.  Values  for  attenuated  long-term 
significant  wave  heights  are  listed  in  Table  24.  These  estimates  represented 
the  best  estimates  for  the  extreme  wave  climate  at  St.  Paul  Harbor,  based  on 
analyses  of  all  the  presently  available  data. 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAOEfl»7l»B  Dim  Bnlmrmd) 


PREFACE 
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was  performed  under  the  direct  supervision  of  Dr.  Dennis  R.  Smith,  Chief,  Pro¬ 
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Mr.  Orson  P.  Smith  and  Ms.  Jane  M.  McKee. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 


(metric)  units  as  follows: 

Multiply 

By 

To  Obtain 

feet 

0.3048 

metres 

miles  (US  nautical) 

1.852 

kilometres 

miles  (US  statute) 

1.609347 

kilometres 

square  miles  (US  statute) 

2589.988 

square  kilometres 

yards 

0.9144 

metres 

8 
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EXTREMAL  ANALYSIS  OF  HINDCAST  AND  MEASURED 
WIND  AND  WAVE  DATA  AT  KODIAK.  ALASKA 


PART  I:  INTRODUCTION 


Purpose  of  Stud; 


I .  this 


1.  This  study  was  conducted  to  provide  an  analysis  of  hindcast  and 
measured  wave  and  wind  data  for  Kodiak,  Alaska.  The  purpose  of  the  analysis 
was  to  obtain  long-term  extreme  wave  conditions  for  the  St.  Paul  Harbor  at 
Kodiak.  This  study  was  accomplished  by  means  of  extremal  analyses  of  deep¬ 
water  hindcast  data  from  the  nearest  Wave  Information  Study  (WIS)  grid  point 
outside  Chiniak  Bay.  Long-term  wind  measurements  were  used  to  validate  the 
results  obtained  from  the  analyses  of  the  hindcast  data.  The  local  wave  cli¬ 
mate  and  attenuation  factor  for  swell  crossing  the  reef  into  St.  Paul  Harbor 
were  derived  using  measured  wave  data  from  the  area.^ 


Physical  Description  of  Kodiak  and  Vicinitv 


2.  The  city  of  Kodiak  is  located  on  the  northeastern  shore  of  Kodiak 
Island,  on  the  western  Gulf  of  Alaska,  about  1,250  air  miles*  northwest  of 
Seattle  and  250  miles  southwest  of  Anchorage.  Kodiak  Island  is  3,588  square 
miles  in  area  and  is  mostly  mountainous  terrain  rising  to  over  4,000  ft  in 
places.  Its  shoreline  is  characterized  by  deep  glacial  fiords  separated  by 
rocky  peninsulas  and  many  smaller  islands.  The  center  of  the  city  lies  on 
the  Kodiak  Island  side  of  a  narrow  channel  defined  by  Near  Island.  The  Port 
of  Kodiak's  deep-draft  facilities  are  southwest  of  the  city  on  the  northwest 
shore  of  that  50-  to  60-ft-deep  area  of  Chiniak  Bay  known  as  St.  Paul  Harbor. 
St.  Paul  Harbor  is  defined  by  a  series  of  small  islands  and  submerged  rocky 
reefs  a  few  feet  deep  extending  from  the  offshore  side  of  Near  Island  2  miles 
to  the  southwest  to  Just  beyond  Puffin  Island.  Further  to  the  southwest  of 
St.  Paul  Harbor  is  Womens  Bay,  the  site  of  the  US  Coast  Guard  Kodiak  Air 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


Station.  Chiniak  Bay,  offshore  of  the  city  and  St.  Paul  Harbor,  is  defined 
by  Cape  Chiniak  and  Long  Island  and  is  exposed  to  the  northern  half  of  the 
Gulf  of  Alaska.  Cape  Chiniak  offers  protection  from  the  Pacific  Ocean  to  the 
south. 

3.  The  specific  area  of  interest  to  this  study  is  the  deep-draft  termi¬ 
nal  operated  by  the  Port  of  Kodiak  on  St.  Paul  Harbor.  This  container  dock  is 
fully  exposed  to  St.  Paul  Harbor,  and  its  operations  are  intermittently  dis¬ 
rupted  by  long-period  swell  which  passes  during  bad  weather  over  the  reefs 
defining  St.  Paul  Harbor. 

4.  Developments  by  the  State  of  Alaska  in  Dog  Bay  on  the  southwest 
side  of  Near  Island  are  under  way;  they  were  proposed  in  1976  by  the  Corps  of 
Engineers  and  are  also  of  interest  in  this  study.  Dog  Bay  is  sheltered  from 
Chiniak  Bay  by  the  southern  tip  of  Near  Island  but  is  exposed  to  occasional 
strong  winds  out  of  Womens  Bay,  which  generate  seas  that  are  hazardous  to 
small  craft.  The  map  of  the  area  given  in  Figure  1,  taken  from  the  Alaska 
Coastal  Data  Collection  Program  Data  (ACDCP)  report,  shows  the  locations  of 
measurement  devices  and  other  features  discussed  in  this  report. 


Figure  1 .  Kodiak  area  map  and  data  collection 
and  telemetry  system 
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PART  II:  HINDCAST  EXTREMAL  ANALYSIS 


Hindcast  Data  Selection 


5.  The  hindcast  wave  data  set  used  in  this  study  was  provided  on  mag¬ 
netic  tape  by  WIS.  The  data  consisted  of  one  record  of  wave  climate  charac¬ 
teristics  including  significant  wave  height,  period,  and  direction  for  every 
3  hr  for  the  years  from  1956  to  1975  for  a  total  of  approximately  58,400  rec¬ 
ords  (Ragsdale  1983).  A  computer  program  was  developed  to  search  this  tape 
for  records  having  specific  directions  and  magnitudes.  The  program  writes 
tabular  listings  of  the  selected  data  records  for  further  data  reduction. 

The  WIS  grid  point  (Point  17)  that  is  nearest  to  the  Chiniak  Bay  area  is 
120  nautical  miles  east  of  Kodiak  (see  Figure  2).  Events  producing  signif¬ 
icant  wave  heights  of  6  m  or  more  were  chosen  and  maximum  significant  wave 
height  and  associated  period  were  then  selected  for  each  event.  Appendix  A 
provides  a  listing  of  the  78  resulting  maximum  significant  wave  heights, 
each  with  its  date  of  occurrence  and  period.  The  choice  of  6  ra  as  a  selec¬ 
tion  criterion  was  arbitrary.  However,  the  choice  was  modified  and  vali¬ 
dated  in  terms  of  the  resulting  extremal  analysis.  The  data  were  surveyed 
for  significant  storm  events  with  wave  directions  corresponding  to  a  direc¬ 
tion  window  defined  by  Cape  Chiniak  to  the  south  and  Long  Island  to  the  north. 
This  window  was  taken  to  be  approximately  95  to  140  degrees  relative  to  True 
North  for  waves  traveling  from  the  north  and  90  degrees  for  waves  from  the 
east.  The  direction  window  was  constructed  using  the  approximate  center  of 
the  St.  Paul  Harbor  area  (a  little  to  the  northwest  of  Puffin  Island)  as  the 
vertex  of  a  triangle  with  the  other  two  vertices  defined  by  Cape  Chiniak  and 
Long  Island. 


Probability  Distribution  for  the 
Number  of  Storms  per  Year 

6.  The  number  of  storms  per  year  producing  significant  wave  heights  of 
6  m  or  more  is  listed  in  Table  1.  Note  that  for  1957  there  was  no  event  pro¬ 
ducing  significant  wave  heights  of  6  m  or  more;  whereas  1963,  1968,  and  1969 
each  produced  seven  such  events.  The  number  of  storms  per  year  is  listed  with 
its  observed  frequencies  in  Table  2.  A  common  assumption  for  studies  of  this 
type  is  that  the  number  of  storms  per  time  interval  is  distributed  according 
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Table  1 


Annual  Number  of  Storms  Producing  Significant 


Wave  Heights  of  6  m  or  More 


Number  of 
Storms 


78  =  N 


Table  2 


Frequency  Table  for  the  Number  of  Storms  per 


1-Year  Interval  Producing  Significant 


Wave  Heights  of  6  m  or  More 


Number  of 
Storms 


Frequency 


to  a  Poisson  probability  distribution  (Borgman  and  Resio  1982).  The  Poisson 
distribution  has  probability  density 


L*  <  •  . 

FS 


uxe~u  n  ,  ,  , 

xi  i  x  ■  0  i  1  (  2  ,  3  i  ... 


The  variable  x  is  the  number  of  storms  per  year  and  the  parameter  u  is  the 
average  number  of  storms  per  year  in  this  case,  and  u  =  78  storms/20  years  or 
u  =  3.9  . 

7.  The  chi  square  goodness  of  fit  test  is  used  to  test  the  validity  of 


the  Poisson  assumption  (Miller  and  Freund  1977).  The  test  is  based  on  the 
2 

statistic  xn 


where 


x  *  z 

n  i=i 


£  \°i  “  Ei> 


0^  =  observed  frequency  of  years  with  i  -  1  storms 

=  Poisson  expected  frequency  of  years  with  i  -  1  storms 


If  xjj  is  smail  when  compared  to  theoretical  chi  square  values,  then  the 


probability  that  the  Poisson  assumption  is  valid  will  be  great.  Table  3 

2 

lists  the  values  for  Oj  and  Et  and  the  resulting  value  for  xn 


Table  3 

Chi  Square  Test  for  the Poisson  Model 


Number  of 
Storms 

°i 

Ei 

(°i  -  Ei)2' 

0 

1 

1  3 

2  5 

0.525 

2 

3 

6  7 

l:So  7-08 

0.001 

4 

5 

l  8 

3.90  c 

3.04 

0.127 

6 

7 

1  - 

^ ’98  3  08 

1.10  ** 

0.275 

X*  =  0.928 

Some  of  the  cells  have  been  combined  in  the  test  to  minimize  the  impact  of 
small  cell  counts  on  the  resulting  statistic.  The  theoretical  chi  square 
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values  are  found  In  The  Handbook  of  Mathematical  Functions  (Abramowitz  and 
Stegun  1972). 


For  a  chi  square  with  5-degrees  of  freedom  (n  =  4) 

Pr^  *  0.71l)  =  0.95 
Pr(xj;  *  1.064)  =  0.90 
and  by  means  of  linear  interpolation 

Pr(x4  *  0.928)  =  0.92 

Since  the  computed  statistic  is  small  compared  to  92  percent  of  all  possible 
chi  square  values,  it  is  concluded  that  the  Poisson  assumption  is  valid.  This 
result  will  be  instrumental  in  defining  such  quantities  as  return  period  and 
nonencounter  probability  in  the  following  sections  of  this  report. 


Extremal  Theory 

8.  There  are  several  theoretical  probability  distributions  that  have  been 
successfully  used  in  the  fitting  and  subsequent  extrapolation  of  extreme  wave 
conditions.  These  are  the  Extremal  or  Fisher-Tippett  Type  I  distribution,  the 
Lognormal  distribution,  the  Log  Extremal  distribution,  and  the  Weibull 
distribution.  These  distributions  have  cumulative  probability  functions  as 
listed  in  Table  4. 

Table  4 

Extremal  Models 


Extremal  Type  I : 


F(x)  -  exp 


|.exp 


-  »  <  x  <  “ 


(3) 


(Continued) 


Lognormal : 


I 


m 


1  f  1  r  1/lnh- 

oh  exP  |_”  2  (  a~~ 


0  <  x  < 


Log  Extremal: 


F(x)  =  exp 


[■  (■?)! 


0  <  x  < 


Welbull: 


F(x)  =  1.0  -  exp 


[-M 


9.  The  theoretical  cumulative  probability  function  is  fit  to  data  by 
means  of  the  plotting  position  formula.  If  the  data  sample  given  by  x^  , 
x2  xn  is  ranked  in  ascending  order  denoted  by  <  Y^2)  <•••<  Y^ 

where  Y^kj  is  called  the  kth  order  statistic,  then  the  plotting  position 
formula 

K  =  rb  <7> 

represents  the  estimate  of  the  data  cumulative  probability  function.  If  this 
is  set  equal  to  the  proposed  theoretical  cumulative  probability  function  F(x) 
from  Table  4  then 

K-rh-  FKk>  *  B1  <8> 

where  A  and  B  are  scale  and  location  parameters,  respectively.  The  in¬ 
verse  of  the  function  in  Equation  8  is 


. . 


(9) 
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If  the  plot  of  F“Uk/n  +  1)  with  Y^)  approximates  a  straight  line  with 
slope  A  and  intercept  B  then  the  proposed  theoretical  distribution  is  ac¬ 
cepted  Sometimes  more  than  one  of  the  possible  distributions  will  yield  a 
stra  ght  line  fit.  In  this  case,  the  better  of  these  is  usually  that  which 

A 

best  fits  the  upper  tail  of  the  function  F^  .  However,  some  subjective 
judgment  is  required  in  such  cases. 

10.  The  quantity  known  as  the  return  period,  R  ,  is  defined  to  be  the 
mean  value  of  the  random  number  of  observations  preceding  and  including  the 
first  exceeds  e  of  a  specified  wave  threshold  x  .  In  terms  of  the  cumula¬ 
tive  probability  function  and  the  Poisson  model  parameter  u 


R  = 


1 


u[1  -  F(x) ] 


(10) 


Another  isefui  measure  is  the  nonencounter  probability,  NE(x)  ,  or  the  prob¬ 
ability  that  for  a  design  life  L  the  largest  wave  condition  is  less  than  or 
equal  to  x  in  value.  For  the  Poisson  model 


NE(x)  =  exp 


(11) 


Note  that  if  L  =  R  then 


NE(x)  =  0.37 


Thus  the  probability  of  encountering  a  condition  larger  than  the  R  year  re¬ 
turn  period  condition  in  R  years  is  0.63.  This  demonstrates  the  misleading 
nature  of  the  return  period  in  that  during  R  years,  there  is  a  63  percent 
chance  of  encountering  an  R  year  extreme  condition.  Care  should  be  taken 
when  using  return  period  as  the  only  criterion  for  extremal  prediction. 


Extremal  Plotting 


11.  A  computer  software  package  called  EXPLOT  was  developed  at  the  US 
Army  Engineer  Waterways  Experiment  Station  (WES)  on  the  Honeywell  DPS1  system 
to  plot  the  formula  given  in  Equation  9.  The  plotting  was  performed  for  each 
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of  the  functions  listed  in  Table  4.  The  Weibull  distribution  was  computed  for 
values  of  C  =  1.0  and  C  =  2.0  from  Equation  6.  These  values  represent  the 
general  range  of  Weibull  shape  parameters  that  are  applicable  to  wave  condi¬ 
tions.  Note  that  if  C  =  1.0  ,  the  Weibull  reduces  to  an  exponential  distri¬ 
bution  and  if  C  =  2.0  ,  the  Rayleigh  distribution  (Petrauskas  and  Aagaard 
1971).  The  slope  and  intercept  from  Equation  9  are  estimated  by  computing  the 
least  squares  linear  regression  line  corresponding  to  Equation  9  (Issacson  and 
Mackenzie  1981).  For  a  listing  of  estimated  values  of  A  and  B  along  with 
the  corresponding  parameter  estimates  for  each  specific  function  see  Table  5: 

Table  5 

Estimated  Parameters  for  the  Hypothesized  Extremal  Models 
_ Model _  A  B  u  _ o 


Extremal  Type 

I 

-6.804 

0.981 

6.936 

1.019 

Lognormal 

-8.058 

1.128 

7.144 

0.886 

Log  Extremal 

-14.952 

7.742 

7.742 

6.898 

Weibull,  C  = 

1.0 

-4.615 

0.745 

6.195 

1.342 

Weibull,  C  = 

2.0 

-1.839 

0.363 

5.066 

2.755 

12.  Appendix  B  contains  the  resulting  data  plots  for  each  of  the  pro¬ 
posed  extremal  models.  The  horizontal  axis  entitled  Cumulative  Probability 
Scale  denotes  the  actual  values  of  the  function  F(x)  from  Table  4  that  cor¬ 
respond  to  the  data  values  on  the  vertical  axis.  The  other  horizontal  axis  is 
the  return  period  as  defined  in  Equation  10.  By  inspection  it  is  seen  that 
the  Lognormal  model  does  not  fit.  The  Weibull  with  C  =  1 .0  ,  or  exponential, 
also  displays  significant  curvature.  The  Log  Extremal  model  shows  less  curva¬ 
ture,  but  the  plot  still  deviates  from  linearity  in  the  upper  and  lower  tails 
of  the  distribution.  The  Extremal  Type  I  and  the  Weibull  with  C  =  2.0 
(Rayleigh)  both  look  fairly  linear  except  for  the  lowest  16  points.  Since 
the  choice  of  6  m  for  the  data  selection  threshold  was  arbitrary,  it  is  in¬ 
tuitively  appealing  to  recompute  the  analysis  without  the  lower  16  points. 

This  results  in  a  total  of  62  storm  events  with  maxima  greater  than  or  equal 
to  6.4  m  and  a  revised  Poisson  parameter  of  u  =  3.1  storms  per  year.  The 
chi  square  goodness  of  fit  test  for  the  Poisson  distribution  of  storms  per 
year  was  recomputed  with  the  reduced  data  set.  The  chi  square  value  was 
=  1.35  and  ^Pr  xjj  *  1.35^  *  0.85  .  The  chi  square  statistic  is  still 
smaller  than  85  percent  of  all  possible  chi  square  values;  therefore,  it  is 
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concluded  that  the  Poisson  model  still  holds.  The  extremal  plots  for  the  re¬ 
duced  data  sets  are  displayed  in  Figures  3  and  4.  The  Extremal  Type  I  and  the 
Weibull  C  =  2.0  both  appear  to  fit  the  largest  62  data  points  well.  The  Ex¬ 
tremal  Type  I  is  the  preferred  model  because  (a)  it  appears  to  fit  better  in 
the  upper  tail  of  the  data,  (b)  it  is  one  of  the  three  possible  extremal  as¬ 
ymptotes  (Borgman  and  Resio  1982),  of  which  the  Weibull  is  not  a  member  and, 
therefore,  has  a  better  theoretical  basis  than  does  the  Weibull,  and  (c)  the 
Extremal  Type  I  is  a  two-parameter  model  whereas  the  Weibull  has  three  param¬ 
eters.  The  Weibull  shape  parameter  C  makes  it  possible  to  fit  the  distribu¬ 
tion  to  data  with  a  higher  degree  of  accuracy  than  for  two-parameter  models. 
This  fact  does  not  necessarily  mean  that  extrapolations  beyond  the  data  extent 
will  be  improved  by  the  higher  precision  fit.  The  extrapolations  could  be 
very  unreliable  if  the  Weibull  were  used  when  a  simpler  model  was  more  appro¬ 
priate.  Results  will  be  given  for  both  models,  but  it  is  believed  that  the 
Extremal  Type  I  is  qualitatively  the  better  model.  Table  6  contains  the  val¬ 
ues  of  A  and  B  and  the  resulting  distribution  parameters.  The  return 
periods  and  associated  significant  wave  heights  were  computed  for  both  models 
and  ar.e  presented  in  Table  7. 
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cumulative  probability  scale 

Figure  3.  Extremal  plot  for  62  largest  maximum 
significant  wave  heights,  Extremal  Type  I 
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Figure  4.  Extremal  plot  for  62  largest  maximum 
significant  wave  heights,  Weibull,  C  =  2.0 


Table  6 

Estimated  Parameters  for  the  Selected  Extremal  Models 


Model 

A 

B 

p 

a 

Extremal  Type  I 
Weibull,  C  =  2.0 

-7.567 

-2.130 

1.036 

0.384 

7.304 

5.547 

0.965 

2.604 

Table  7 

Predicted  Significant  Wave  Heights  and  Associated 
Return  Periods  for  the  Selected  Models 


Return  Period 
years 


5 

10 

20 

40 


Extremal  Type  I 


Significant 
Wave  Height,  m 


9.92 

10.60 

11.28 

11.95 
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Table  7  (Concluded) 


Return  Period 
years 


Significant 
Wave  Height,  m 


Extremal  Type  I  (Continued) 


50 

60 

100 


5 

10 

20 

40 


12.17 

12.35 

12.84 


Weibull,  C  =  2.0 

9.86 

10.37 

10.84 

11.26 


50 

60 

100 


11.39 

11.50 

11.78 


13.  The  concept  that  we  refer  to  as  "return  period"  has  been  shown  to 
be  prone  to  misinterpretation  (Borgman  1963).  The  statement  "50-year  wave" 
may  be  taken  to  mean  that  waves  of  magnitude  greater  than  or  equal  to  the 
50-year  value  occur  on  the  average  of  once  every  50  years.  The  phrase  "on  the 
average"  is  the  key  to  interpreting  this  concept.  As  was  seen  in  the  discus¬ 
sion  of  nonencounter  probabilities,  the  chance  of  encountering  a  50-year  wave 
in  any  given  50-year  interval  is  about  63  percent.  The  nonencounter  probabil¬ 
ities  given  in  Table  8  are  defined  to  be  the  probabilities  for  any  given  de¬ 
sign  life  L  that  an  R  year  wave  is  not  encountered.  These  values  can  be 
used  to  give  further  meaning  to  the  return  periods  and  associated  significant 
wave  heights  reported  in  Table  6. 


Table  8 

Estimated  Nonencounter  Probabilities 
for  the  Extremal  Type  I  Model 


Design  Life, 

Return 

Period, 

years  (R) 

years  (L) 

20 

40 

50 

60 

100 

20 

0.37 

0.61 

0.67 

0.72 

0.82 

40 

0.14 

0.37 

0.45 

0.51 

0.67 

50 

0.08 

0.29 

0.37 

0.43 

0.61 

60 

0.05 

0.22 

0.30 

0.37 

0.55 

100 

0.01 

0.08 

0.14 

0.19 

0.37 
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Estimate  Reliability 


14.  Methods  for  computing  the  reliability  of  extreme  wave  condition 
estimates  can  be  found  in  various  reports  and  articles  (Petrauskas  and  Aagaard 
1971,  Isaacson  and  Mackenzie  1981,  Borgman  1983).  The  method  given  by  Borgman 
requires  relatively  little  computation  while  the  other  methods  either  do  not 
apply  to  prediction  or  entail  extensive  numerical  simulations.  The  method 
used  here  has  been  successfully  applied  to  the  extremal  analysis  of  hindcast 
data  on  the  coast  of  California  (Borgman  1982). 

15.  The  first  step  in  computing  the  estimate  reliability  is  to  calcu¬ 
late  an  upper  bound  on  the  error  (standard  deviation)  aR  from  estimating 
the  function  F(x)  from  Table  4.  This  upper  bound  is  approximated  by  the 
Kolmogorov-Smirnov  (KS)  limit  for  the  quantity 


|Fr  -  F(x) 


(12) 


A 

for  Fr  from  Equation  7.  The  KS  upper  bound  is  the  value  D(a,n)  such 
that 


PR 


Jmax)FR  -  F(x) )  >  D(a,n)J  =  a 


for  a  =  a  small  number.  In  other  words,  the  number  D(a,n)  is  such  that  the 
maximum  absolute  difference  between  the  cumulative  probability  function  F(x) 

A 

and  its  estimated  value  FR  will  be  greater  than  D(a,n)  only  with  small 
probability  a 
error  where 


The  value  D(a,n)  is  then  assumed  to  represent  a  *  ctf 


Za/2  is  the  standard  normal  variate  such  that  Pr(Z  <  Za/2) 


a/2  ;  thus,  an  approximate  value  for  the  error 


cjp  is  computed  to  be 


-  p(a»n) 


"a/2 


(13) 


For  samples  of  n  >  35  and  a  =  0.05  an  approximate  value  for  D(a,n) 
given  n  =  62  is 


D(0.05,n)  = 


^2 


and 
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Z0.025  =  1,96 


The  value  for  op  represents  the  approximate  sampling  error  or  standard  de¬ 
viation  involved  in  the  estimation  of  F(x)  .  In  order  to  state  this  error 
in  terms  of  the  Y-year  significant  wave  height,  it  is  necessary  to  define 
the  quantity  known  as  the  coefficient  of  variation  or  CV  .  The  CV  is 
the  standard  deviation  of  a  variable  divided  by  the  mean  of  the  variable. 

At  the  median  where  F(x)  =  0.5  ,  the  CV  for  the  quantity  F(x)  will  be 
approximately 


CVp  =  53 


The  CV  for  F(x)  is  approximately  equal  to  the  CV  for  Fy(x)  »  the  cumu¬ 
lative  probability  density  for  the  Y-year  significant  wave  height  (Borgman 

1982);  thus,  if  cu  is  the  standard  deviation  for  F„(x)  and  since  Fp(x) 
FR  y  « 

=0.37  for  the  R-year  wave,  then 

aP  =  0.37CV- 


=  0.07 


16.  The  final  step  in  the  reliability  analysis  is  to  express  the  error 
in  terms  of  the  R-year  significant  wave  height  error  o  .  This  is  accom- 

XR 

plished  by  means  of  the  approximate  formula 


where  aD  =  the  slope  of  FR  evaluated  at  XF 


The  formula  for  aR  is 


3F(Xr) 


exp  |uR[ 1  -  F(Xr)]| 


. 


for  u  =  the  Poisson  parameter  or  the  mean  number  of  storm  events  per  year. 
The  slope  aR  for  the  Extremal  Type  I  is  approximately 


3R  =  ,03 


(20) 


for  all  values  of  R  from  20  to  100. 

17.  The  error  due  to  the  hindcast  model  can  also  be  included  in  the 
reliability  analysis.  The  WIS  hindcast  model  when  compared  to  measured  data 
shows  at  most  a  difference  of  about  0.04  for  percent  occurrence  of  signifi¬ 
cant  wave  heights  larger  than  6  m  (Figure  11,  Corson  and  Resio  1981).  If 
this  value  is  assumed  to  be  the  error  in  estimating  F(x)  due  to  the  hind¬ 
cast  model  then  <jp  from  Equation  14  becomes 


op  =  0.09  +  0.04  =  0.13 


(21) 


Recomputing  Equations  15  through  17,  the  result  is 


oP  =  0.10  (22) 

FR 

and  then  by  Equation  18  for  the  Extremal  Type  I 


oY  =  0. 10 

*R 


(23) 


18.  The  error  in  the  extremal  prediction  to  a  very  rough  approximation 
is  about  one-tenth  of  a  meter.  It  is  not  yet  known  how  accurate  this  error 
estimation  technique  is,  and  its  results  should  be  used  with  some  subjective 
consideration. 


Wave  Period  Distribution 

19.  The  Joint  probability  distribution  for  significant  wave  height  ver¬ 
sus  period  is  listed  in  Table  9.  For  the  extreme  wave  condition,  the  periods 
tend  to  be  larger  than  12  sec  with  periods  at  14.3  sec  for  the  three  most  ex¬ 
treme  conditions  in  the  20-year  period.  This  tendency  suggests  the  use  of  a 
period  of  near  14.3  sec  for  the  computation  of  design  conditions. 
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Table  9 

Joint  Frequency  Distribution  of  Wave  Period 
Versus  Significant  Wave  Height 


Significant 
Wave  Height 
_ m _ 

6  <  Hs  <  7 

7  <  Hs  <  8 

8  <  Hs  <  9 

9  <  Hs  <  10 

10  <  Hs  <  11 

11  <  Hs  <  12 


11- 

6 

4 

0 

0 

0 

0 


11-12 

14 
9 
4 
0 
0 
0 


Wave  Period,  sec 


12-13 

9 
2 
8 
5 
1 
0 


13-14 
0 


14-15 

3 
5 

4 
1 
2 
1 


20 


m 

V 

i 


mz&mzxm immm. 


I  IN  P*J T "J  1*. 


*.1  ■  m.in.1  Ul  K  J  i"i«  l  *  U'i.i  uu'Air.^  u'.i  wn  i!  *  if'rrs.'ip.^uR  mmj  ■  v  ^  «—*■-*  'u<w«  u*»>«  «•  «  i 


PART  III:  ANALYSIS  OF  MEASURED  WAVE  DATA 

Data  Description 

20.  Data  were  provided  to  WES  by  US  Array  Engineer  District,  Alaska 
(NPA),  on  magnetic  tape  duplicating  that  published  in  Data  Reports  1  and  2  of 
the  ACDCP.  These  data  consisted  of  three  sets  of  tables:  "Summary  Table  for 
Wind/Wave  Data,"  "Summary  Table  for  Wind  Data,"  and  Energy  Spectrum"  tables 
presenting  data  collected  by  two  Waver ider  buoys.  The  Waverider  buoys,  de¬ 
noted  as  buoy  1  (outer  buoy)  and  buoy  2  (inner  buoy),  are  located  in  Chiniak 
Bay  and  St.  Paul  Harbor,  and  are  in  77  m  and  16  m  of  water,  respectively  (see 
Figure  1).  The  data  provided  included  measurements  taken  from  October  1981 
through  September  1983.  Also  provided  on  magnetic  tape  was  an  edited  version 
of  National  Weather  Service  (NWS)  anemometer  data  from  1945  through  1982  mea¬ 
sured  by  an  instrument  located  between  the  Kodiak  Airport  and  Womens  Bay  at 
the  present  site  of  the  US  Coast  Guard  Kodiak  Air  Station.  The  detailed  for¬ 
mat  of  the  original  data  and  modifications  made  during  the  analysis  will  be 
described  in  the  following  paragraphs. 

Data  Selection 


21.  Extreme  wave  conditions  for  the  north  Pacific  tend  to  occur  dur¬ 
ing  the  winter  months.  This  is  exemplified  by  the  frequency  distribution 
for  months  with  extreme  wave  conditions  producing  significant  wave  heights 
of  6  m  or  more  given  in  Table  10.  The  values  in  Table  10  were  computed  using 
the  hindcast  data. 

Table  10 

Frequency  Distribution  of  Months  During  Which  Extreme 
Events  Occurred  for  the  20-Year  Hindcast  Period 


Month 

Frequency 

October 

9 

November 

17 

December 

20 

January 

22 

February 

7 

March 

2 

April 

1 

21 


1 


22.  The  winter  "storm  season"  is  defined  to  be  the  months  from  October 
through  April.  With  this  definition,  it  is  possible  to  consider  only  those 
months  that  would  be  expected  to  produce  extremes  in  the  measured  data.  This 
eliminates  possible  biases  due  to  the  mixture  of  monthly  maxima  resulting  from 
storms  with  those  resulting  from  regular  conditions. 


Extremal  Plotting  of  Monthly  Maxima 


23.  The  monthly  maxima  were  fit  to  the  Weibull  probability  model  with 
shape  parameter  C  =  1  ,  and  C  =  2  and  to  the  Extremal  Type  I.  As  for  the 
hindcast  data,  the  Extremal  Type  I  displays  a  better  fit  than  does  the 
Weibull.  Figures  5  through  10  represent  the  extremal  plots  for  the  measured 
data.  The  extreme  wave  conditions  and  associated  return  periods  given  in 
Table  11  were  computed  from  the  fitted  Extremal  Type  I  model.  The  values 


Table  11 


Significant  Wave  Height  and  Return 


Period  for  Measured  Data 


Return  Period 


Significant 
Wave  Height 


Inner  Buo\ 


Outer  Buoy 
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CUMULATIVE  PROBABILITY  SCALE 
HE I BULL »  C  r  1.0 

Figure  10.  Extremal  plot  for  monthly  maximum  signifi¬ 
cant  wave  heights  measured  at  the  Outer  Buoy  location, 

Weibull  C  =  1.0 


given  in  Table  11  are  intended  for  comparison  with  the  hindcast  values  but  not 
for  use  in  constructing  long-term  extreme  wave  conditions.  It  is  a  generally 
accepted  rule  of  thumb  that  any  extrapolation  beyond  twice  the  extent  of  the 
data  becomes  highly  unreliable  (Borgman  and  Resio  1982);  hence,  2  years  of 
data  may  be  useful  for  predicting  2  to  5  years  of  wave  conditions  but  are  not 
very  reliable  for  5-  to  100-year  predictions. 

Discussion 


24.  The  use  of  measured  data  to  make  long-term  wave  condition  estimates 
has  been  discussed  in  several  recent  articles  (Carter  and  Challenor  1978, 
Challenor  1982,  and  Wang  and  LeMehaute  1983).  The  key  aspect  of  the  theory  of 
extremes  is  that  in  looking  at  enough  years  of  data,  the  chance  of  seeing  most 
types  of  hazardous  events  is  good.  In  this  respect,  measured  data  will  become 
useful  as  longer  records  are  made  available.  At  this  time,  hindcast  data  ful¬ 
fill  the  need.  The  wave  data  that  were  measured  at  Kodiak  provide  an  inter¬ 
esting  comparison  to  hindcast  results. 

25.  It  is  surprising  that  the  outer  buoy,  or  buoy  1,  did  not  display  at 
least  one  maximum  significant  wave  height  in  exceedance  of  6  m.  If  it  is  as¬ 
sumed  that  the  hindcast  grid  point  represents  the  same  general  wave  condition 
as  the  outer  buoy,  and  since  there  were  78  events  producing  significant  wave 
heights  of  6  m  or  more  in  the  20-year  hindcast  time  period,  the  probability 

of  seeing  two  consecutive  years  with  no  significant  wave  heights  greater  than 
6  m  at  the  outer  buoy  is  nearly  zero.  Because  of  this,  it  is  not  possible  to 
either  support  or  contradict  the  results  from  the  hindcast  analysis  using  the 
measured  wave  data.  If  consecutive  hindcast  and  measured  values  were  avail¬ 
able,  this  behavior  could  be  studied;  however,  such  is  not  the  case  at  this 


PART  IV:  ENERGY  ATTENUATION  FACTOR 


Analytical  Procedures 


26.  Two  FORTRAN  computer  programs  provided  the  principal  tools  for 
analysis  of  the  data  measured  at  Kodiak:  "Statistical  Package  for  the  Social 
Sciences"  (SPSS)  and  a  program  titled  "WINSUM."  SPSS  is  a  comprehensive  soft¬ 
ware  system  of  advanced  statistical  routines,  originally  developed  by  the 
National  Opinion  Research  Center  of  the  University  of  Chicago.  SPSS  statisti¬ 
cal  routines  are  highly  versatile  and  applicable  to  most  scientific  and  engi¬ 
neering  data.  WINSUM  was  developed  by  the  Coastal  Engineering  Research  Center 
for  analysis  of  periodic  wind  measurements.  These  two  programs  and  a  number 
of  supportive  editing  and  plotting  routines  were  executed  on  the  WES  Honeywell 
DPS-1  system. 

27.  Data  from  the  "Summary  Tables  for  Wind/Wave  Data"  and  the  "Summary 
Tables  for  Wind  Data"  were  first  reduced  to  a  simple  listing  of  consecutive 
individual  readings  in  a  predetermined  FORTRAN  format.  Individual  runs  of 
SPSS  accomplished  additional  grouping  value  corrections,  accounting  of  missing 
values,  and  labeling  required  by  particular  statistical  procedures.  Zero  val¬ 
ues  were  included  with  missing  values.  Similarly,  the  NWS  wind  data  were  re¬ 
duced  to  a  simple  listing  of  consecutive  readings  in  a  format  compatible  with 
WINSUM  and  correlations  related  to  the  physical  characteristics  of  the  mea¬ 
surement  site  were  also  accomplished. 

Wave  Data  Analysis 

28.  The  "Summary  Table  for  Wind/Wave  Data"  presents  the  significant 
(zero-moment)  wave  heights  and  peak  periods  for  each  sample  from  both  wave 
buoys.  A  description  of  how  these  values  were  computed  from  the  original  time 
series  data  is  provided  by  the  ACDCP  (1983a,  b).  The  complete  record  was  ana¬ 
lyzed  for  the  frequency  distribution  of  wave  heights  in  classes  of  2  ft  each 
and  of  peak  periods  in  classes  of  2  sec  each.  The  results  of  this  procedure 
are  presented  in  Tables  12-15. 

29.  Dependence  of  conditions  measured  at  the  inner  buoy  on  conditions 
measured  concurrently  at  the  outer  buoy  were  first  tested  by  SPSS  scattergram 
routines  which  included  simple  linear  regression.  These  measurements  were 


Table  12 

Frequency  Distribution  Analysis  for  the  Outer  Buo' 
Location  Significant  Wave  Height  Measurements 


Frequency 


Category 

Label 

Code 

Absolute 

Relative 

percent 

Adjusted 

percent 

Cumulative 

percent 

Less  than  2  ft 

1 

1,889 

36.3 

52.8 

52.8 

2  to  4  ft 

2 

1,054 

20.3 

29.5 

82.3 

4  to  6  ft 

3 

378 

7.3 

10.6 

92.8 

6  to  8  ft 

4 

147 

2.8 

4.1 

97.0 

8  to  10  ft 

5 

60 

1.2 

1.7 

98.6 

10  to  12  ft 

6 

29 

0.6 

0.8 

99.4 

12  to  14  ft 

7 

17 

0.3 

0.5 

99.9 

14  to  16  ft 

8 

1 

0.0 

0.0 

99.9 

Over  16  ft 

9 

2 

0.0 

0.1 

100.0 

-1 

1,625 

31.2 

Missing 

100.0 

Total 

5,202 

100.0 

100.0 

Mean 

1.772 

Standard  error 

0.018 

Mode 

1.000 

Standard  deviation 

1.088 

Kurtosis 

5.108 

Skewness 

1.978 

Minimum 

1.000 

Maximum 

9.000 

Median  1 .447 
Variance  1.183 


Valid  cases  3,577 
Missing  cases  1,625 


Table  13 

Frequency  Distribution  Analysis  for  the  Inner  Buoy 
Location  Significant  Wave  Height  Measurements 


Frequency 


29 


Category 

Label  Code 

Absolute 

Relative 

percent 

Adjusted 

percent 

Cumulative  j 

percent 

Less  than  2  ft 

1 

4,062 

78.1 

83.1 

83.1 

2  to  4  ft 

2 

765 

14.7 

15.6 

98.7  P 

* 

4  to  6  ft 

3 

62 

1.2 

1.3 

100.0  ■' 

6  to  8  ft 

4 

2 

0.0 

0.0 

100.0  , 

-1 

311 

6.0 

Missing 

100.0 

Total  5,202 

100.0 

100.0 

! 

! 

• 

* 

* 

j 

Mean  1.183 

Standard  error 

0.006 

Median 

1 . 102  1 

Mode  1 .000 

Standard  deviation  0.421 

Variance 

0.177  | 

J 

Kurtosis  4.459 

Skewness 

2.215 

Range 

3.000 

Minimum  1.000 

Maximum 

4.000 

1 

Valid  cases  4,891 

1 

V 

• 

• 

s 

I 

• 

j 

Missing  cases  311 

! 

v 

> 

i 

S 

a 


Table  14 


Frequency  Distribution  Analysis  for  the  Outer  Buoy 


Location  Peak  Period  Measurements 


Category 


Label 

Code 

Absolute 

Less  than  2  sec 

1 

14 

2  to  4  sec 

2 

373 

4  to  6  sec 

3 

454 

6  to  8  sec 

4 

736 

8  to  10  sec 

5 

1,161 

10  to  12  sec 

6 

379 

12  to  14  sec 

7 

210 

Over  14  sec 

8 

250 

-1 

1,625 

Total 

5,202 

Relative 

percent 

Adjusted 

percent 

Cumulative 

percent 

0.3 

0.4 

0.4 

7.2 

10.4 

10.8 

8.7 

12.7 

23.5 

14.1 

20.6 

44.1 

22.3 

32.5 

76.5 

7.3 

10.6 

87.1 

4.0 

5.9 

93.0 

4.8 

7.0 

100.0 

31.2 

Missing 

100.0 

100.0 

100.0 

y  •  v 
>» 

Mean 

4.645 

Standard  error 

0.027 

Median 

4.682 

& 

Mode 

5.000 

Standard  deviation 

1.600 

Variance 

2:559 

:K 

Kurtosis 

-0.284 

Skewness 

0.231 

Range 

7.000 

* 

Minimum 

1.000 

Maximum 

8.000 

r-  -* 
i  »  . 


Valid  cases  3,577 
Missing  cases  1,625 


Table  15 

Frequency  Distribution  Analysis  for  the  Inner  Buo 


Location  Peak  Period  Measurements 


Category 
_ Label 

Less  than  2  sec 
2  to  4  sec 
4  to  6  sec 
6  to  8  sec 


Code  Absolute 


1  79 

2  615 

3  782 

4  1,112 


_ Frequency _ 

Relative  Adjusted 


ercent 


ercent 


Valid  cases  4,891 

Missing  cases  311 


Cumulative 
ercent 


8  to  10  ; 

sec 

5 

1,223 

23.5 

25.0 

77.9 

10  to  12 

sec 

6 

404 

7.8 

8.3 

86.2 

12  to  14 

sec 

7 

206 

4.0 

4.2 

90.4 

Over  14  : 

sec 

8 

470 

9.0 

9.6 

100.0 

-1 

311 

6.0 

Missing 

100.0 

Total  5,202 

100.0 

100.0 

Mean 

4.466 

Standard  error 

0.025 

Median 

4.3' 

Mode 

5.000 

Standard  deviation  1.757 

Variance 

3.01 

Kurtosis 

-0.386 

Skewness 

0.401 

Range 

7.01 

Minimum 

1.000 

Maximum 

8.000 

taken  for  approximately  20  min  at  the  outer  buoy  immediately  followed  by  a 
similar  sample  of  the  inner  buoy.  This  sequence  and  the  time  lapse  involved 
is  effectively  the  same  as  concurrent  measurement,  since  conditions  at  sea, 
particularly  more  extreme  conditions,  very  rarely  change  that  fast.  The  scat- 
tergram  of  HS2  (inner  buoy  significant  wave  height)  and  HS1  (outer  buoy)  shows 
some  linearity,  as  indicated  in  Figure  11  (R  =  0.39).  TP2  (inner  buoy  peak 
period)  and  TP1  (outer  buoy)  do  not  show  any  statistical  relation,  though  vi¬ 
sually  the  scattergram  (Figure  12)  indicates  TP2  tends  to  equal  TP1  in  many 
cases.  Attempts  to  run  scattergrams  of  HS22  with  HS12  and  HS2/gTP22  with 

p 

HSI/gTPI  produced  no  meaningful  information.  The  significant  wave  heights 
were  originally  computed  as  representative  of  energy  in  the  sea  state,  so 
the  relation  of  energy  levels  is  inherent  in  any  comparison  of  HS2  with  HS1. 

p 

Squaring  these  values  apparently  only  exaggerated  the  scatter.  H/gTc,  as  a 
measure  of  wave  steepness,  showed  no  statistical  trends  apparently  because 
the  peak  periods  by  themselves  didn’t  either. 

30.  Attempts  to  relate  measured  wave  conditions  to  concurrently  mea¬ 
sured  wind  direction  did  not  yield  any  dramatic  conclusions,  but  they  did  re¬ 
inforce  the  intuitive  assumption  that  the  most  severe  wave  conditions  exist 
with  winds  from  Chiniak  Bay.  Figure  13  is  a  scattergram  of  HS2  with  DIR  (di¬ 
rection),  in  which  cases  with  HS1  <  1.0  m  and  TP1  <  4.0  sec  were  excluded. 
The  grouping  of  lower  HS2  values  in  the  range  280-360  deg  indicates  the  rela¬ 
tive  frequency  of  offshore  winds.  HS2  values  above  1  m  are  grouped  in  the 
range  45-180  deg.  Figure  14  is  a  scattergram  of  HS1  with  DIR,  showing  simi¬ 
lar  indications.  The  ratio  TP2/TP1  plotted  against  DIR  in  Figure  15  shows 
the  tendency  for  TP2  to  equal  TP1  with  directional  grouping  similar  to  the 
wave  height  scattergrams.  Womens  Bay  winds  blow  toward  Puffin  Island  in  the 
200-  to  250-deg  range.  This  range  has  few  points  plotted  on  the  three  scat¬ 
tergrams,  indicating  severe  wave  conditions  concurrent  with  winds  from  Womens 
Bay  are  rare. 

31.  SPSS  procedures  attempting  to  correlate  wind  velocity  at  Puffin 
Island  with  significant  wave  heights  did  not  reveal  any  meaningful  conclu¬ 
sions.  The  relationship  of  wind  velocity  to  significant  wave  height  is  known 
to  be  highly  nonlinear,  however,  and  includes  a  number  of  other  equally  crit¬ 
ical  independent  variables.  Wind  velocity  data  are,  therefore,  treated  sep¬ 
arately  in  this  report  with  a  view  toward  application  of  these  nonlinear 
relationships. 
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Figure  14.  Scattergrajn  of  significant  wave  height  at  the 
Outer  Buoy  (HS1)  versus  wind  direction  (DIR) 


32.  The  most  significant  statistical  relationship,  following  prelimi¬ 
nary  analyses,  appeared  to  be  the  reduction  of  significant  wave  height  between 
buoy  1  and  buoy  2,  as  revealed  by  the  scattergrams .  SPSS  procedures,  which 

2 

computed  multiple  stepwise  regression  of  HS2  with  factors  including  HS1,  HS1  , 
HS13,  TP1,  and  TP12  were  executed.  Table  16  is  a  summary  of  a  run  where  cases 
with  TP1  <  4  sec  were  excluded.  The  independent  variables  are  listed  in 
order  of  their  relative  significance  in  the  regression  equation: 

V  .  A  .  B,X,  .  B2X2  .....  BnXn  (24) 

where  A  and  Bn  are  constants  and  Xn  is  the  independent  variable.  The 
changes  in  the  correlation  (Rc)  indicate  that  HS1  provides  the  only  important 
correlation,  with  TP1  as  next  most  important  though  comparatively  insignifi- 

p 

cant.  Nonlinear  correlations  were  not  revealed  by  inclusion  of  the  HS1  , 

HS13,  or  TP12  terms. 


Table  16 

Regression  Summary  for  the  Multiple  Linear  Regression  of  HS2 
with  HS1 .  HS12  HS13  TP1,  and  TP12  for  TP1  >  4  sec 


Multiple 

R 

RSQ 

Simple 

Var iable 

R 

Square 

Change 

R 

B 

Beta 

HS1  -  Sig  Wave  Ht  1 

0.74649 

0.55725 

0.55725 

0.74649 

0.44623 

1 . 1 1870 

TP1  -  Peak  Period  1 

0.76011 

0.56266 

0.00542 

-0.02561 

-0.02978 

-0.30323 

TP1SQ  -  TPeakl  Squared 

0.75218 

0.56578 

0.00312 

-0.05501 

0.00108 

0.24627 

HS1SQ  -  HSIG1  Squared 

0.75358 

0.56788 

0.00209 

0.67612 

-0.07965 

-0.66127 

HS1CU  -  HSIG1  Cubed 

0.75513 

0.57022 

0.00235 

0.56785 

0.01090 

0.33109 

(Constant) 

0.26593 

33.  Scattergrams  of  HS2  with  HS1  which  excluded  cases  with  TP1  values 
below  an  increasing  lower  limit  showed  increasing  correlation.  Figure  16  is 
a  scattergram  of  HS2  with  HS1  where  cases  with  TPA  <  12  sec  were  excluded. 
The  R2  computed  was  0.73  with  405  cases  included.  The  grouping  near  axes 
intersection  and  along  the  HS2  =  HS1  line  is  typical  of  all  the  scatter¬ 
grams  where  only  lower  TP1  values  were  excluded.  Figure  17  is  a  scattergram 
of  HS2  with  HS1  where  the  following  conditions  were  specified:  TP1  >  12  sec  , 
HS2  >  0.5  m  ,  and  HS1  >  1.5  m  .  The  R2  computed  was  0.75  with  38  cases 


included.  The  complete  regression  statistics  for  this  condition  are  presented 
in  Table  17.  The  linear  regression  equation  resulting  from  this  analysis  is: 

HS2  =  0.003  +  0.374  HS1  (25) 

The  intercept  is  notably  very  near  zero  indicating  that  an  attenuation  factor 
of  0.37  for  the  conditions  stated  above  has  significance.  In  fact  the  90  per¬ 
cent  confidence  interval  on  the  intercept,  given  by  (-0.17,  0.18),  contains 
zero  indicating  that  the  intercept  is  not  significantly  different  from  zero 
in  a  statistical  sense.  The  slope  or  attenuation  factor  has  90  percent  con¬ 
fidence  interval  (0.31,  0.43). 

Table  17 

Linear  Regression  Summary  for  Inner  Buoy 
Significant  Wave  Height  (HS2)  and  Outer 
Buoy  Significant  Wave  Height  (HS1) 


Variable(s)  Entered  on  Step  Number  1 


Multiple  R 

R  Square 

Adjusted  R  Square 
Standard  Error 

0.86888 

0.75495 

0.74814 

0.19925 

Analysis  of 
Variance 

DF 

Sum  of 
Squares 

Mean 

Square 

F 

Regression 

Residual 

1 

36 

4.40304 

1.42921 

4.40304 

0.03970 

110.90697 

Variables  in  the  Equation 

Variable 

B 

Beta 

Standard 
Error  B 

F 

HS1 

0.37425 

0.86888 

0.03554 

110.907 

(Constant)  0.00311 


PART  V:  ANALYSIS  OF  MEASURED  WIND  DATA 


NWS  Data 

34.  The  NWS  anemometer  is  located  at  the  Coast  Guard  Station  within 
Womens  Bay.  The  anemometer  is  approximately  100  yd  from  the  water  line  at 
mean  sea  level,  elevation  10  m.  Prior  to  1973,  a  Navy  anemometer  was  at  the 
same  location,  but  at  an  elevation  of  5  m.  The  NWS/Navy  anemometer  provided 
hourly  data  from  November  1945  through  December  1982. 

Data  adjustment 

35.  Using  the  method  presented  in  the  Coastal  Engineering  Technical 
Note  CETN-I-5,  "Method  of  Determining  Adjusted  Windspeed,  Ua  ,  for  Wave 
Forecasting,"  WES  personnel  adjusted  the  data.  The  adjustments  included  the 
following: 

a.  Correction  to  10-m  level  for  data  prior  to  1973. 

b.  Correction  to  hourly  averages  from  5-min  averages  every  hour. 

c.  Adjustment  of  overland  readings  to  overwater  readings. 

d.  Correction  for  nonconstant  coefficient  of  drag. 

e.  Correction  for  instability  due  to  air-sea  temperature  differ¬ 
ences  for  directions  where  the  fetch  is  greater  than  10  miles 
(an  unstable  condition  was  assumed,  since  no  temperature  data 
were  available). 

These  corrections  were  made  so  that  the  wind  data  could  be  applied  directly  to 
wave  forecasting  curves. 

Wind  summary 

36.  A  computer  program  was  developed  to  produce  summaries  of  the  hourly 
winds  and  of  the  3-,  6-,  8-,  and  10-hourly  winds.  The  summaries  relate  wind 
speed  and  duration  to  recurrence  probability  for  16  directional  sectors.  The 
maximum  annual  wind  was  also  found  for  use  in  an  extremal  analysis  to  be  pre¬ 
sented  later. 

37.  The  procedure  to  compute  the  N-hour  average  wind  was  as  follows: 

a.  Delete  bad  observations  (less  than  1  percent  of  the 
325,050  observations  were  bad). 

b.  Compute  arithmetic  average  of  all  wind  speeds  in  consecutive 
N-hour  intervals  to  estimate  wind  speed  (calm  values  were 
counted  as  zero  wind  speeds). 

c.  Compute  vector  sum  of  all  noncalm  observations  in  the  N-hour 
interval. 


I 

j 

d.  Compute  direction  of  vector  sum  from  c  to  estimate  wind 
direction. 

38.  The  distribution  of  the  hourly  winds  with  respect  to  direction  is 
shown  as  a  wind  rose  (Figure  18).  The  wind  rose  shows  the  distribution  of 
winds  in  16  directional  sectors  as  a  percent  of  all  winds.  The  bars  in  each 
direction  are  broken  into  10-mph  intervals. 

39.  The  distributions  of  wind  speed  and  duration  relating  to  return 
period  are  presented  in  Figure  19.  Figure  19  shows  combined  data  from  all 
sectors  (0-360  deg).  The  curves  represent  the  visual  best  fit  to  the  data. 

From  these  curves,  the  wind  speed  for  a  given  return  period  can  be  obtained. 

In  general,  these  curves  should  not  be  extrapolated  beyond  twice  the  period 
of  record  (37  years). 

Puffin  Island  Data 


HO.  The  Puffin  Island  anemometer  is  part  of  the  ACDCP.  Puffin  Island 
is  located  1.5  miles  offshore  in  Chiniak  Bay.  The  anemometer  is  located  at  an 
elevation  of  23  m  (island  elevation)  plus  6.1  m  (elevation  of  sensor  above  the 
island).  The  anemometer  supplied  about  16  months  of  usable  data  starting  late 
in  1981.  The  anemometer  was  not  working  in  October  and  November  of  1981  nor 
in  October  of  1982,  months  when  major  storms  are  expected.  The  anemometer  was 
also  out  during  other  extreme  events,  so  the  data  for  extreme  events  are  un¬ 
reliable.  The  observations  were  1-hr  averages  at  3-hr  intervals. 

Data  adjustments 

Hi.  Again,  the  raw  data  were  adjusted  using  the  method  presented  in 
CETN-I-5.  The  adjustments  included: 

a.  Correction  to  10-m  level. 

b.  Adjustment  of  overland  readings  to  overwater  readings. 

c  Correction  for  nonconstant  coefficient  of  drag. 

d.  Correction  for  instability  due  to  air-sea  temperature  differ¬ 
ences  for  directions  where  the  fetch  is  greater  than  10  miles 
(an  unstable  condition  was  assumed,  since  no  temperature  data 
were  available). 

A  calibration  correction  factor  of  1.06  was  also  applied  as  recommended  by  the 
Alaska  District.  An  additional  adjustment  was  made  to  account  for  the  influ¬ 
ence  of  the  island  height  and  shape  on  the  airflow  (Simiu  and  Scanlon  1978). 


This  correction  is 
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b  a 


where 


in  which 


Lin  (- 
\z  > 


h  =  height  of  the  hill 
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L  =  factor  of  island  length 

zQ  =  roughness  length 

u°  =  dimensionless  quantity  representing  the  perturbation  of 
the  upwind  velocity  due  to  the  presence  of  the  hill 

z  =  anemometer  height  above  the  island 

l  =  thickness  of  the  internal  boundary  layer 

This  correction  assumes  that  L  is  much  larger  than  h  and  that  it  is  rural 
terrain.  A  maximum  correction  of  O.65  was  applied  in  the  direction  of  the 
short  axis  of  the  island,  and  a  minimum  correction  of  0.85  was  applied  in  the 
direction  of  the  long  axis  of  the  island.  These  corrections  were  made  so  wind 
data  could  be  applied  directly  to  wave  forecasting  curves. 

Wind  summary 

42.  The  same  computer  program  used  to  analyze  the  NWS  data  was  used  for 
the  Puffin  Island  data.  Data  were  available  for  only  the  3-hr  averages,  so 
summaries  were  produced  for  the  3-  and  6-hr  winds.  Of  the  4,800  observations, 
17  percent  were  bad. 

43.  The  distribution  of  the  hourly  wind  with  respect  to  direction  is 
shown  as  a  wind  rose  (Figure  20).  The  wind  rose  shows  the  distribution  of 
winds  in  16  directional  sectors  as  a  percent  of  all  winds.  The  bars  are 
broken  into  10-mph  intervals. 

44.  The  short  period  of  record  for  the  data  from  Puffin  Island  anemom¬ 
eter  makes  it  difficult  to  create  reliable  curves  for  wind  speed  and  duration 
versus  return  period.  The  best  use  for  this  data  is  for  comparison  with  the 


15%  WIND  ROSE 


Figure  20.  Wind  rose  plot  for  Puffin  Island  data 


Comparison  of  NWS  and  Puffin  Island  Data 


45.  The  Puffin  Island  data  seem  to  compare  well  in  velocity,  direc¬ 
tion,  duration,  and  return  period  with  the  NWS  data.  Figure  21  shows  the 
distribution  of  wind  speed  and  duration  relative  to  return  period  (3-  and 
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6-hr  duration).  The  solid  curves  represent  the  distribution  for  the  NWS  data 
(see  Figure  19).  The  Puffin  Island  data  follow  the  NWS  distribution  closely 
in  the  lower  parts  of  the  curves.  This  is  the  region  where  the  Puffin  Island 
data  are  most  reliable.  This  implies  that  the  NWS  wind  data  are  a  good  esti¬ 
mate  of  the  Puffin  Island  wind  conditions  which  would  give  a  good  estimate  of 
the  open  water  conditions  in  St.  Paul  Harbor  and  Chiniak  Bay  at  Kodiak.  The 
following  section  presents  the  extremal  analysis  of  the  winds  using  the  NWS 
wind  data.  The  extremal  curves  are  recommended  for  design  purposes. 


Analysis  of  Annual  Maxima  for  NWS  Data 


46.  The  NWS  data  for  1-,  3-,  6-,  8-,  10-,  12-,  15-,  18-,  21-,  24-,  27-, 
and  30-hr  wind  speed  averages  were  searched  for  yearly  maxima  for  input  to  the 
extremal  plotting  routine,  EXPLOT.  The  resulting  plots  are  provided  in  Appen¬ 
dix  C.  It  is  apparent  for  all  but  the  1-hr  averages  that  the  smaller  eight 
data  points  have  been  smoothed  out  and  can  no  longer  be  considered  extremes. 
For  this  reason  the  points  in  question  were  not  used  when  fitting  the  least 
squares  lines  to  the  data,  and  therefore,  do  not  affect  the  resulting  extrap¬ 
olated  wind  speeds.  Table  18  contains  a  full  listing  of  the  predicted  wind 
speeds  for  given  return  periods  and  hourly  averages.  The  analysis  was  com¬ 
puted  for  all  directions,  the  offshore  exposure  (45-180  deg  relative  to  0  deg 
from  North)  and  Womens  Bay  (202.5-247.5  deg)  separately  as  indicated  in  Appen¬ 
dix  C  and  Table  18.  The  Extremal  Type  I  or  Fisher-Tippett  I  distribution  was 
used  in  this  analysis.  The  Fisher-Tippett  or  Extremal  Type  II  distribution 
results  were  also  computed  and  found  to  be  unsatisfactory.  In  general,  the 
chosen  extremal  model  fits  the  data  very  well  as  is  demonstrated  by  the 
plots. 


Wave  Forecasts  Based  on  NWS  Wind  Data 

47.  The  long-term  (extremal)  distributions  of  wind  velocities  averaged 
over  varying  durations  for  all  directions  for  247.5-292.5  deg  (Womens  Bay), 
and  for  45-157.5  deg  (Gulf  of  Alaska)  are  presented  in  Appendix  C.  This  in¬ 
formation  can  be  used  to  derive  equivalent  significant  wave  heights  and  peak 
spectral  periods,  given  appropriate  assumptions  of  water  depth  and  fetch. 

Two  sectors  of  interest  were  defined  which  correspond  to  the  7-mile  fetch  of 


Table  18 

Maximum  Windspeeds,  Miles  per  Hour  with 
Hours  Averaged  and  Return  Period 


Hours 

Return  Period,  years 

Averaged 

5 

10 

20  40  50 

60 

100 

All  Directions 


1 

94.0 

103.4 

112.4 

121.3 

124.1 

126.4 

132.9 

3 

80.1 

89.2 

98.0 

106.5 

109.2 

111.5 

117.7 

6 

74.0 

83.1 

91.8 

100.3 

103.1 

105.3 

111.5 

8 

70.9 

79.4 

87.5 

95.5 

98.0 

100.1 

105.9 

10 

68.4 

76.8 

85.0 

92.9 

95.4 

97.5 

103.3 

12 

49.3 

56.2 

62.3 

68.3 

70.2 

71.8 

76.1 

15 

46.2 

52.3 

58.2 

64.0 

65.9 

67.4 

71.6 

18 

44.8 

50.5 

56.1 

61.5 

63.2 

64.6 

68.6 

21 

42.9 

48.4 

53.6 

58.7 

60.4 

61.7 

65.4 

24 

41.5 

47.2 

52.6 

57.9 

59.6 

61.0 

64.8 

27 

37.8 

41.8 

45.6 

49.3 

50.5 

51.5 

54.2 

30 

36.9 

41.4 

45.7 

50.0 

51.3 

52.4 

55.5 

Offshore  Exposure 

1 

72.2 

79.7 

86.8 

93.8 

96.0 

97.8 

102.9 

3 

62.1 

70.7 

78.9 

87.0 

89.5 

91.6 

97.5 

6 

56.1 

63.5 

70.7 

77.8 

80.0 

81.9 

87.0 

8 

53.0 

59.5 

65.8 

72.0 

74.0 

75.6 

80.1 

10 

50. 1 

56.1 

62.0 

67.6 

69.4 

71.0 

75.0 

Womens  Bay 

1 

52.5 

61.7- 

70.5 

79.1 

81.8 

84.1 

90.4 

3 

37.2 

47.1 

56.7 

66.0 

69.0 

71.5 

78.3 

6 

26.5 

31.1 

35.5 

39.7 

41.1 

42.2 

45.3 

8 

26.1 

31.5 

36.7 

41.8 

43.5 

44.8 

48.5 

10 

24.3 

29.8 

35.1 

40.2 

41.9 

43.2 

47.0 

Womens  Bay  to  Dog  Bay, 

and  the 

1  virtually  unlimited 

fetch  across 

i  the  Gulf 

of 

Alaska. 

The  depths  along  the 

Womens  Bay 

fetch  vary 

from  a  few 

feet  to  over 

50  ft. 

Shallow  water 

forecasting  curves 

(Vincent  and  Lockhart 

1983)  for 

an 

assumed  constant  35-ft  depth  and  7-mile  fetch  were  used  to  translate  veloc¬ 
ities  and  durations  at  selected  return  periods  from  Figures  C1-C5  to  equiva¬ 
lent  significant  wave  heights  and  peak  periods  for  the  Womens  Bay  fetch,  as 
presented  in  Table  i9.  The  fetch  limitation  makes  the  higher  velocities  of 
1-hr  duration  the  extreme  condition  in  terms  of  wave  height. 
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Table  19 

Wave  Forecast  for  Womens  Bav  Fetch  Based  on  NWS  Wind  Data 


Duration 

hr 


_ Return  Period.  ; 

10  20  40  50 


1.8,  4.2  2.0,  4.2  2.2,  4.7  2.6,  5.0  2.6,  9.0  2.7,  5.0  2.9,  5.2 

1.1,  3.7  1.4,  4.0  1.7,  4.3  2.1,  4.6  2.3,  4.7  2.3,  4.8  2.5,  4.9 


Note:  First  number  under  Return  Period  indicates  Significant  Wave  Height,  m; 
second  number  indicates  Peak  Period,  sec. 

48.  A  similar  translation  of  winds  from  the  Gulf  of  Alaska  (Figures  C6- 
C17)  was  performed  using  the  deepwater  forecasting  curves  with  duration  as  the 
only  limitation.  Table  20  indicates  a  duration  of  24  hr  as  the  extreme  condi¬ 
tions  in  this  case  with  low  probability  deepwater  significant  wave  heights 

Table  20 

Wave  Forecast  Tor  Gulf  of  Alaska  Based  on  NWS  Wind  Data 


Duration 

hr 


12 

5-5, 

10.9 

6.4, 

11.4 

7.3, 

12.0 

8.2, 

12.6 

8.5, 

12.8 

8.8, 

13.0 

9.1, 

13.3 

15 

5.8, 

11.7 

6.7, 

12.3 

7.9, 

13.1 

8.8, 

13.8 

9.1, 

14.0 

9.5, 

14.1 

10.4, 

14.8 

18 

6.7, 

12.6 

7.6, 

13.4 

8-5, 

14.1 

9.5. 

14.7 

10.1, 

15.0 

10.4, 

15.1 

11.0, 

15.5 

21 

7.0, 

13.3 

8.2, 

14.2 

9.1, 

15.0 

10.4, 

15.8 

10.7, 

15.9 

11.0, 

16.1 

11.9, 

16.6 

24 

7.0, 

13.9 

8.5, 

14.6 

10.1, 

15.8 

11.3, 

16.4 

11.6, 

16.8 

11.9, 

16.9 

12.8, 

17.3 

27 

7.0, 

14.2 

8.2, 

14.8 

8.8, 

15.5 

9.8, 

16.1 

10. 1, 

16.3 

10.4, 

16.4 

11.0, 

17.0 

30 

6.7, 

13.8 

8.2, 

15.3 

9.8, 

16.2 

11.0, 

17.0 

11.3, 

17.1 

11.6, 

17.2 

12.5, 

17.9 

Note:  First  number  under  Return  Period  indicates  Significant  Wave  Height,  m;  second  number 
indicates  Peak  Period,  sec. 


very  close  to  those  predicted  using  the  WIS  data  base  (see  Table  21).  The 
duration  over  21  hr  for  the  low  probability  events  have  corresponding  fetches 
on  the  order  of  300  to  400  nautical  miles.  Though  no  physical  barriers  to  the 
wind  exist,  this  distance  is  on  the  scale  of  the  synoptic  weather  patterns. 

The  variations  of  velocity  across  weather  patterns  of  this  scale  are  modeled 
by  the  WIS  data  base,  but  are  not  using  the  NWS  point  source  of  data.  This 
inadequacy,  along  with  the  other  measurement  errors  inherent  in  land-based 
anemometers,  indicates  that  the  statistics  from  the  WIS  data  base  are  more 
reliable,  as  applied  here,  for  wave  conditions  at  Kodiak. 
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Table  21 

Comparison  of  Extremes  Predicted 
from  WIS  Hindcast  Data  and 
NWS  Wind  Data  Analysis 


Return  Period 

Significant  Wave 
Height,  m 

years 

WIS 

NWS 

5 

9.92 

7.0 

10 

10.6 

8.5 

20 

1 1 .28 

10.1 

40 

11.95 

11.3 

50 

12.17 

11.6 

60 

12.35 

11.9 

100 

12.84 

12.8 

52 


... 


PART  VI:  MATER  DEPTH  AND  HAVE  BREAKING  CONSIDERATIONS 


49.  The  outer  buoy  location  water  depth  is  about  77  m.  For  linear 
waves  with  periods  of  12  sec  and  longer,  this  does  not  constitute  deep  water. 
The  deepwater  assumption  holds  for  h/Lo  >0.5  where  h  =  water  depth  and 
Lo  =  deepwater  wavelength.  Values  for  h/Lo  with  h  =  77  m  and  periods 
of  12  sec  and  longer  are  h/Lo  <0.34  ;  thus,  the  waves  traveling  from  deep 
water  will  be  somewhat  affected  by  the  water  depth.  Corrections  for  this 
are  available  in  wave  tables  (Skovgaard  et  al.  1974)  that  report  the  ratio 
of  wave  height  to  deepwater  wave  height  as  a  function  of  h/Lo  .  Values  for 
this  ratio  are  listed  in  Table  22  for  the  range  of  wave  periods  of  interest. 
Wave  breaking  limits  were  computed  using  a  controlling  depth  of  7.5  m  for  the 
area  between  Puffin  Island  and  Popov  Island.  It  was  found  that  for  the  ex¬ 
treme  waves  of  this  study  the  attenuation  factor  does  not  need  to  be  modified 
using  an  upper  limit  of  significant  wave  height. 


Wave  Period 

H 

sec 

Ho 

14.0 

0.932 

14.3 

0.929 

14.5 

0.929 

14.8 

0.926 

15.0 

0.923 

15.3 

0.920 

15.5 

0.920 

15.8 

0.918 

16.0 

0.916 

16.3 

0.916 

16.5 

0.914 

16.8 

0.914 

17.0 

0.913 

50.  The  peak  wave  period  that  most  accurately  represents  the  local 
weather  conditions  comes  from  the  WIS  hindcast  data  set.  The  WIS  hindcast 
period  for  the  20-year  event  is  14.3  sec.  This  value  is  also  sufficient  for 


the  50-  and  100-year  wave  conditions.  The  ratio  of  wave  height  in  77  m  of 
water  to  deepwater  wave  height  for  waves  with  14. 3-sec  periods  is  0.929; 
thus,  the  50-  and  100-year  deepwater  significant  wave  heights  given  in 
Table  7,  12.17  and  12.84,  respectively,  will  be  11.3  and  11.9  m  at  the 
outer  buoy  location. 


PART  VII:  SUMMARY  AND  CONCLUSIONS 


51.  The  WIS  hindcast  data  applicable  to  Chiniak  Bay  were  surveyed 
for  extreme  wave  conditions.  The  resulting  62  significant  wave  heights  were 
found  to  fit  the  Extremal  Type  I  distribution.  The  resulting  extrapolated 
significant  wave  heights  are  reported  in  Table  7.  The  associated  reliability 
analysis  indicates  that  if  the  Extremal  Type  I  is  the  proper  model,  then  the 
standard  error  in  estimating  the  extremes  is  about  one-tenth  of  a  meter.  It 
should  be  stressed  that  the  major  source  of  error  is  due  to  the  choice  of  the 
model  and  is  thus  not  included  in  this  number. 

52.  Measured  data  from  the  area  were  also  found  to  fit  an  Extremal 
Type  I  distribution.  The  results  from  this  analysis  tend  to  underestimate 
long-term  conditions  such  as  the  10-  to  50-year  events.  This  is  mostly  due 
to  the  short  data  record  length  (2  years).  Extremes  computed  using  the  NWS 
wind  data  are  similar  to  those  given  by  the  hindcast  analysis.  For  the 
50-year  event  they  differ  by  about  six-tenths  of  a  meter,  and  for  the 
100-year  event  they  are  approximately  the  same. 

53. -  The  peak  wave  period  that  most  accurately  represents  the  local 
weather  patterns  is  that  computed  by  the  WIS  hindcast.  The  wave  period  asso¬ 
ciated  with  the  20-year  extreme  event  is  14.3  sec.  This  value  should  be  suf¬ 
ficient  for  the  50-  and  100-year  events  as  well.  The  deepwater  to  intermedi¬ 
ate  water  depth  (77  m  at  the  outer  buoy)  wave  height  correction  for  14. 3-sec 
periods  is  0.929.  This  correction  was  applied  to  the  significant  wave 
heights  listed  in  Table  7  and  the  values  listed  in  Table  23  were  obtained. 

Table  23 

Depth  Corrected  to  Outer  Buoy  (77  m)  Significant 
Wave  Heights  Versus  Return  Period 


Return  Period 

Significant  Wave 

years 

Height,  m 

5 

9.22 

10 

9.85 

20 

10.48 

40 

11.10 

50 

11.31 

60 

11.47 

100 

11.93 

55 


54.  The  deep-  to  shallow-water  attenuation  factor  from  Equation  25  is 
0.374.  The  values  in  Table  23  multiplied  by  this  factor  result  in  the  inner- 
buoy  or  shallow  water  extreme  values  given  in  Table  24. 


Table  24 


Attenuated  Significant  Wave  Heights 


Versus  Return  Period 


Return  Period 


Significant  Wave 
Height,  m 


5 

10 

20 

40 


3.45 

3.68 

3.92 

4.15 


50 

60 

100 


4.23 

4.29 

4.46 


The  values  given  in  Table  24  along  with  the  wave  period  of  14.3  sec  are  the 
estimates  for  the  long-term  extreme  events  pertaining  to  the  St.  Paul  Harbor 
area,  based  on  analyses  of  all  the  currently  available  data. 
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Kodiak  Storm  Data 


Period 


Date 

HSMAX.  m 

sec 

01-18-56 

6.20 

12.50 

01-28-38 

8.80 

12.50 

02-12-58 

6.60 

11.10 

02-24-58 

6.90 

11.10 

11-20-58 

7.80 

11.10 

12-10-58 

6.30 

10.00 

01-16-59 

11.70 

14.30 

01-21-59 

7.20 

10.00 

11-16-59 

7.40 

11.10 

01-31-60 

9.90 

14.30 

02-07-60 

8.90 

14.30 

10-24-60 

7.50 

11.10 

11-22-60 

7.00 

11.10 

12-20-60 

6.70 

11.10 

01-19-61 

9.20 

12.50 

01-24-61 

6.20 

11.10 

11-26-61 

6.30 

11.10 

10-27-62 

8.10 

12.50 

11-07-62 

6.30 

11.10 

11-23-62 

7.20 

10.00 

11-28-62 

6.30 

10.00 

12-30-62 

6.00 

11.10 

02-08-63 

8.40 

12.50 

10-21-63 

6.80 

11.10 

10-31-63 

9.30 

12.50 

11-11-63 

6.70 

14.30 

11-22-63 

6.50 

12.50 

11-25-63 

7.20 

11.10 

12-21-63 

8.50 

11.10 

01-12-64 

6.90 

11.10 

01-24-64 

6.60 

10.00 

01-20-64 

9.40 

12.50 

12-01-64 

8.20 

12.50 

10-02-65 

6.30 

11.10 

10-05-65 

7.60 

11.10 

01-13-66 

7.30 

12.50 

01-26-66 

8.60 

12.50 

11-01-66 

7.40 

11.10 

01-10-67 

7.10 

14.30 

02-27-67 

6.00 

12.50 

(Continued) 


A2 


5JI 

hi 


12-02-67 

12-07-67 

12-23-67 

01-16-68 

01-22-68 


03-21-68 

11- 04-68 
11-18-68 
11-28-68 

12- 14-68 


02-07-69 

02-11-69 

10- 11-69 

11- 19-69 

12- 02-69 


12-19-69 

12-23-69 

01-20-70 

10- 31-70 

11- 11-70 


01-17-71 

01-15-72 

03-09-72 

11-25-72 

12-17-72 


12-24-72 

01-27-73 

12-06-73 

12-13-73 

12-19-73 


12-23-73 

12-27-73 

01-17-74 

04-13-74 

10-30-74 


01-12-75 

11- 19-75 

12- 21-75 


Kodiak  Storm  Data  (Concluded) 


8.80 

9.10 

8.00 

6.30 

9.10 

6.60 

6.70 

7.20 

10.20 

7.00 


10.10 

7.80 

6.10 

6.30 

8.60 


7.10 

10.00 

8.00 

6.10 

8.40 


Period 


12.50 

12.50 

14.30 

12.50 

10.00 


14.30 

11.10 

11.10 

10.00 

12.50 


10.00 

11.10 

10.00 

11.10 

12.50 


12.50 

12.50 

11.10 

11.10 

12.50 


14.30 

12.50 

11.10 

14.30 

14.30 

14.30 

14.30 

12.50 

12.50 

14.30 


10.00 

12.50 

14.30 

11.10 

14.30 


14.30 

12.50 

11.10 


RETURN  PER I 00  IN  TERRS 
10  20  40  60 


90.0  95.0  98.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
LOCNORHAL 


RETURN  PER  100  (N  YEARS 
5  10  20 


80.0  90.0  95.0  98.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
HE I BULL .  C  =  1.0 


.4 


APPENDIX  C 


:  EXTREMAL  DISTRIBUTIONS  OF  WIND  VELOCITIES 
FOR  WOMENS  BAY  AND  GULF  OF  ALASKA 


RETURN  PER I 00  IN  YEARS 
5  10  20  40  60  100 


I  50.0  80.0  90.0  95.0  9B.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
EXTREMAL  TYPE  I 

Figure  Cl.  Womens  Bay,  1-hr  averages 


RETURN  PER I 00  IN  YEARS 
5  tO  20  40  60 


u  20.0 
s 


50.0  SO.O  90.0  95.0  98.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
EXTREMAL  TYPE  I 

Figure  C2.  Womens  Bay,  3-hr  averages 


50.0  80.0  90.0  95.0  98.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
EXTREMAL  TYPE  I 

Figure  C3.  Womens  Bay,  6-hr  averages 


RETURN  PER  100  IN  YEARS 
5  10  20  40  60  100 


o  30.0 

IftJ 


>  10.0 


5 

2  0.0 


SO.O  80.0  90.0  95.0  98.0  99.0 

cumulative  probability  scale 

EXTREMAL  TYPE  I 


Figure  C4.  Womens  Bay,  8-hr  averages 
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RETURN  PER 100  IN  YEARS 


1.0  SO.O  BO .0  90.0  95. 0  9B.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
EXTREMAL  TYPE  I 

Figure  C9.  Gulf  of  Alaska,  8-hr  averages 


RETURN  PERIOO  IN  YEARS 


99.9 


TWENTY  ONE  HR  RVERRCE  HRX  HIND  SPEED  IN  HPH  EIGHTEEN  HR  RVERROE  HRX  WIND  SPEED  IN  HPH 


80.0 


RETURN  PER I 00  IN  TERRS 
5  10  20  40  60  100 


)  SO.O  80.0  90.0  95.0  98.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
EXTREMAL  TYPE  I 

Figure  C13.  Gulf  of  Alaska,  18-hr  averges 


RETURN  PERIOD  IN  YEARS 
S  10  20  40  60  100 


)  SO.O  .  80.0  90.0  95.0  98.0  99.0 

CUMULATIVE  PROBABILITY  SCALE 
EXTREMAL  TYPE  I 

Figure  C14.  Gulf  of  Alaska,  21-hr  averages 


RETURN  PER I 00  IN  YEARS 
S  10  20  40  60  100 


50.0  60.0  90.0  95.0  90.0  99.0 

CUMULATIVE  PR00ABILITY  SCALE 
EXTREMAL  TYPE  I 


Figure  Cl 7.  Gulf  of  Alaska,  30-hr  averages 
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Figure  CM.  Womens  Bay,  8-hr  averages 
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